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Abstract — For networked systems, the control law is typically 
subject to network flaws such as delays and packet dropouts. 
Hence, the time in between updates of the control law varies 
unexpectedly. Here, we present a stability theorem for nonlinear 
model predictive control with varying control horizon in a 
continuous time setting without stabilizing terminal constraints 
or costs. It turns out that stability can be concluded under the 
same conditions as for a (short) fixed control horizon. 

I. Introduction 

In recent years, networked control systems (NCS) received 
growing popularity due to their lower implementation costs 
and greater interoperability compared to standard control 
systems, cf. [8]. NCS designs have been implemented in 
different areas such as robotics, automotive and aeronautical 
applications, see, e.g., [11], [12]. On the backside, however, 
the stability and performance analysis of feedbacks designed 
for networked control systems is more complex, cf. [18]. 

In this paper we investigate stability and performance of a 
nonlinear model predictive controller (MPC) in a prediction 
consistent network scheme. By now, MPC has been under- 
stood quite well even in the nonlinear case, see, e.g. [4], 
[15]. The beauty of the method lies in its simplicity: First, 
an optimal control is computed over a finite optimization 
horizon. Then, a portion of this control is implemented — 
throughout this paper the length of this portion is called 
control horizon — and last the optimization horizon is shifted 
in time rendering the method iteratively applicable. Predic- 
tion consistency formalizes the equivalence of the control 
input history for the actuator and the controller. Due to 
computing and transmission times, delays and dropouts, this 
equivalence is subject to different time levels within the NCS 
components. Hence, the control histories do not coincide 
automatically but this property can be forced using certain 
communication schemes, cf. [1], [6], [14]. At the center of 
these schemes is the requirement of long control sequences 
being transmitted from the controller to the actuator. MPC 
is ideally suited to obtain such a sequence since it not only 
delivers such a control by construction, but its internal model 
can also be used to compensate for time delays between the 
NCS components. 

The prediction consistent approach naturally leads to the 
fact that the control horizon varies over time. The goal 
of this paper is to provide the theoretical foundation for 
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ensuring stability of MPC without stabilizing terminal costs 
or constraints for such varying control horizons in a nonlinear 
continuous time setting. To this end, we extend results from 
[16] where a stability condition was introduced for a fixed, 
typically short control horizon to the networked context. In 
particular, the condition allows us to prove stability for a 
large range of possibly time varying control horizons in 
analogy to [7] where the same assertion was shown for 
discrete time systems. This main result is complemented by 
an analysis of structural properties with respect to the over- 
shoot bound and the decay rate of the imposed controllability 
condition. 

The paper is organized as follows. In Section [II] the 
detailed problem formulation is stated. The results from [16] 



are briefly summarized in Section III and then extended to the 



time varying case in Section IV Finally, we show numerical 



results to illustrate our results and draw some conclusions. 

II. Problem setting 

We consider a nonlinear control systems governed by a 
differential equation 



x(t) = f(x(t),u{t)). 



(1) 



Here, x(t) and u(t) denote the state and the control at 
time t > 0, respectively. The system dynamics is given 
by / : R™ x W n — > W l , state and control constraints are 
represented by suitable subsets X C R™ and U C R m , 
respectively. The trajectory which emanates from initial state 
xq and is manipulated by the control function u : K>o — ► K m 
is denoted by x(t;xo,u). We call a control function u 
admissible for x on the interval [0, T) if the conditions 

x{t;x,u) € X,t£ [0,T], and u(t) G U,t € [0,T), (2) 

hold which is denoted by u £ U x ([0,T)). Furthermore, a 
control function u : R>o — > K m is said to be admissible for 
x if, for each T > 0, u e U x ([0,T)) holds for its restriction 
to [0,T). Then, we write u € Z4([0,oo)). 

Here, we consider a networked situation, that is System ([T| 
is connected to an external controller via a network which 
may be subject to delays and packet dropouts. Within the 
network, we suppose the clocks at sensor, controller and 
actuator to be synchronized, see [19, Section III.C] for a 
relaxation of this assumption. These clocks are denoted by 
t s , t c and t a , respectively. 

Our goal is to stabilize System ([T} at an equilibrium x* for 
which we suppose that a control input u* exists satisfying 
f(x*,u*) = 0. To solve this task in an optimal fashion, we 



introduce continuous running costs I 
which satisfy 



^>o 



l(x*,u*) = and inf l(x,u) > rj(\\x — x*\\),x ^ x* , 

for a /Coo -function V- As usual, a continuous function r/ : 
M> — » K>o is said to belong to /Coo if Tj(0) — 0, it is strictly 
increasing, and unbounded. Utilizing the running costs I, we 
define the cost functional 



<^oo (*£> ^v 



l(x(t; x, u), u(t)) dt 



which we wish to minimize for a given initial value i = ioe 
K n . The corresponding optimal value function is denoted by 



Ko(x) := inf Joo(x,u). 

«GW x ([0,oo)) 



(3) 



Since optimal control problems on an infinite time horizon 
are computationally hard, we use model predictive control 
(MPC) in order to approximately solve this task. To this 
end, we firstly set x := xq. Then, our MPC scheme consists 
of the following three steps: 

• Compute a minimizing control function u* for the 
optimal control problem on a truncated and, thus, finite 
prediction horizon T depending on x, i.e. 

f T 
min Jt(x,u) = min / l(x(t;x,u),u(t))dt. 

ueWi[0,T) u&Uit[0,T)jQ 

(4) 
. Define the MPC feedback law fi T : [0, T) x E™ -*■ W n 
by nr{t,x) := u*(i), t e [0,T). Then, for given 
control horizon 6 G [0, T], implement the first piece 
u*(t)| t6 [o.5) of the computed control at the plant in 
order to obtain 

x^ TS {8;x) = x(S;x,u*). 

Whenever we want to reflect the control horizon 5 in 
the notation we write /j,x,s = Ht- 

• Shift the prediction (optimization) horizon forward in 
time by 5 and obtain the new state measurement x which 
coincides in nominal MPC with x^ T S (S; x). 

Iterative application of this procedure generates a solution 
on the infinite time horizon. The resulting input signal 
and trajectory at time t are denoted by /i MVC (t;xo) and 
x MPC (t;xo), respectively. Here, we tacitly assume that Prob- 
lem Q is solvable and the minimum is attained in each step 
of the proposed MPC algorithm. For a detailed discussion of 
feasibility issues we refer to [4], [15]. 

In order to compensate for the mentioned network flaws, 
we utilize a communication protocol which transmits time 
stamped information between the network components: For 
one, the sensor transmits the latest state measurement x to- 
gether with the measurement time instant t s to the controller. 
Secondly, the controller sends predicted control functions 
Mt.A.ct = m*[o,a) with 5 < A <T where creM corresponds 
to the time at which the control should be applied. Here, 
we assume that the controller computes the input functions 
Mt,A,o- at a predefined sampling rate (typically with sampling 
time « A) based on the most recent measurement available 



in the buffer. To store these controls, we add a buffer to the 
actuator and to the controller, cf. Fig. [T] 
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Fig. 1 . Scheme of the networked control system 

The idea of introducing time stamps is the following: 
Given the latest measurement time instant t s , we can use the 
synchronized clocks to compute the transmission delay t sc 
from sensor to controller. Now, in order to be implementable, 
we require that a in the computation of Ht,a,<t is chosen 
such that the control function Ht.a.o arrives at the actuator 
buffer at a time t a < a. To accomplish this, we need to 
know the computing and transmission delays r sc , r c and 
r ca . Since the latter delays are not known at computation 
time, bounds r c max and r ca max are imposed which gives 
us a = t s + t sc + r c max + r ca max , see also Fig. [2] for a 
schematical sketch. Now, we use a model based predictor 
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Fig. 2. Timeline connections within the scheme 

within our controller as in [1], [13], [14], [17] to compensate 



for these delays. If r c + r ca > t c 



holds, the 



packet is considered to be lost. 

While the actuator buffer allows us to store /iT.A.cr in order 
to compute the control input 



u(t) = Ai T ^A lC r(i - o\ x(a)) 



(5) 



the role of the controller buffer is to establish the following 
consistency property: 

Definition 1: (i) A feedback control ^t.a.<t(', %{&)) is 
consistently predicted if the control U[t s .cr) used for the 
prediction of x(a) in the controller is identical to the control 
u [t 3 .a) applied by the actuator. 

(ii) A networked control scheme is prediction consistent if 
at each time f g R the applied control in the actuator ([5]) 
is well defined, i.e. t < a + A, and ht.a.o- is consistently 
predicted. 

The concept of prediction consistency allows us to separate 
the analysis of the feedback law, e.g., in terms of stability 
and performance, from the influence of the network on the 
closed-loop, e.g., in terms of robustness, cf. [2]. Examples 



for prediction consistent network arcitectures can be found, 
e.g., in [1], [6], [14], 

Note that MPC is ideally suited for such an application. 
For one, a model of the plant is already at hand, i.e. the 
prediction of the state measurement x can be done by 
evaluating the system dynamics ([TJ using the control stored 
in the controller buffer. Secondly, the outcome of the MPC 
algorithm is already an optimal control defined on a control 
horizon of length T. Hence, for A < T a control signal of 
length A is readily available for transmission. 

The delays r sc and r ca and thus the length 5 of the portion 
of each fir, A, a used at the actuator may vary with time 
depending on the current network load. In principle, one 
could make this value independent of time by setting 8 = A, 
however, for robustness reasons it is desirable to always use 
the most recent control input, i.e., to keep 8 as small as 
possible [2]. 

Hence, time varying control horizons 5 should be con- 
sidered and our goal is thus to find a condition which 
ensures asymptotic stability of the MPC closed-loop in this 
situation. To this end, we aim at employing the value function 
Vr(-) as a common Lyapunov function and show the relaxed 
Lyapunov inequality 



V T (x^ T (S;x)) <Vt(x) - a / l(x^ T (t]x)^ T {t,x))dt 

Jo 

(6) 

with a £ (0,1] for each feasible state x £ X and all 8 £ 

[t c ,T — t c ], cf. [5] for a discrete setting. 

III. Recap 

Our main tool in order to establish the relaxed Lyapunov 
Inequality (|6]> is a stability condition introduced in [16] for 
fixed control horizon 8. Assumption [2] is needed to state the 
respective result. 

Assumption 2: Let C > 1 and /i > be given. Suppose 
that, for each x £ X and t £ M>o, a control function u x £ 
^[0,00) exists which satisfies 

l{x Ux {t;x),u x (t)) < Ce^*min/(a;,u) =: Ce~^l*(x). 

Assumption [2] is an exponential controllability condition 
in terms of the stage cost with overshoot C and decay rate 
p. Then the main result deduced in [16] reads as follows: 

Theorem 3: Suppose that 77 ( 1 1 .-z; — x*\\) < l*(x), x £ X, 
with r\ £ /Cqo an d Assumption 13] hold. Furthermore, let T > 
8 > and a £ (0, 1) be chosen such that ar,s given by 



1 



0^ _ 1) 



{etiTS)_iy 



( C MT _ !) c _ ^8 _l)V ( eM T - 1) C_ ( e MT-<5) _ 1) c 

(7) 
satisfies the condition q>t,s > cL Then, for each xo £ X, 
the MPC closed-loop solution xt^ s (-;Xo) is asymptotically 
stable and satisfies the suboptimality bound 



Inequality |8]l gives a performance estimates which com- 
pares the resulting MPC closed-loop costs with the theoret- 
ically achievable minimal costs on the infinite time horizon. 
Here, the mono tonicity of the optimal value function Vr(-) 
in the prediction horizon T — an inherent property of 
unconstrained MPC schemes — is crucial in order to deduce 
this bound on the, in general, unknown quantity V^o(-). Note 
that, for given suboptimality index a £ (0, 1), the stability 
condition <xt,s > <5 always holds for a sufficiently large 
prediction horizon T, cf. [16, Section 4.1]. 

IV. Results 

As pointed out in the previous section, in the networked 
context we would like to have a stability criterion for time 
varying control horizon 8. Unfortunately, the stability con- 



■Tin,) = / l(x™%x ),p™ c (t;x ))dt < ^o)_ 



(8) 



dition presented in the previous Section III assumes a fixed 
control horizon. In this section we show how to extend the 
result to varying 8. In particular, we show that the inequality 
&t,s > 5 for some control horizon 8 £ (0, T) implies 
ckT,(5 > 5 for every 5 £ [mm{5, T — 5}, max{<5, T— 5}], This 
result is the key to derive a stability theorem for time varying 
control horizons without imposing additional assumptions. 

To this end, we show symmetry and monotonicity prop- 
erties of the performance bound q>t,s with respect to the 
control horizon 6. Indeed, symmetry is a direct consequence 
of the presented formula. 

Corollary 4: The performance index a-r.s given by For- 
mula |7| satisfies cct,s = a T,T-8, i- e - &T.5 is symmetric with 
symmetry axis 6 = T/2. 

In contrast to the symmetry property given in Corollary 
H deducing monotonicity of Formula (|7j in 6 on (0, T/2] 
is not straightforward. To this end, we apply results from 
[5] and [20]. In the first of these references a discrete 
time counterpart of Theorem [3] was presented based on a 
corresponding version of Assumption[2] In the second article, 
the relation between both the discrete and continuous time 
approaches was investigated. To this end, so called iterative 
refinements were used, i.e. a sequence of discretizations such 
that each element is a partition of its predecessor. Then, 
monotone convergence of the corresponding suboptimality 
bounds to ax.s given by Formula |7]i was shown for the 
discretization parameter tending to zero. We combine these 
results to show that monotonicity carries over from the 
discrete to the continuous time setting. 

Proposition 5: For 5 £ (0, T/2), the performance es- 
timate ot.T.8 given by Formula |7]i has the monotonicity 
property 

a T ,s<a T5 for ~5 E [6, T/2}. (9) 

Proof: We show |9| for T,S,5 £ Q. The assertion for 
real values T,S,S £ M. then follows from the fact that Q is 
dense in R and that ar,s is continuous in T and 5. 

Consider T,8,6 £ Q with 6 £ [8, T/2}. Since 5 and 6 
are rational numbers, there exist r £ Q and m,fh,N £ N 
such that itlt — 8, fhr — 8 Nt — T hold, i.e. r is a 
common denominator. Then, using the abbreviations Nk := 



2 N, to/;, := 2 k m, we define a^ = afe(iV, m) by 



N k N k 

n (7? - 1) n (7.f - 1) 



i- 



nifc + l 



A'fe-nifc+l 



iw-nw 

rn k + l m k +l 



1) 



117? 

N k -m k + l 



N, 



n(7f-i 

7V fc -m fc +l 

(10) 
with 7s fc := C^^cr 2 " with a := e -" T G (0,1). 
Analogously, a^ = ak(N,fh) is defined with m^ := 2 k fh 
instead of mfc. Summarizing, otk and a.k only deviate in the 
parameters m and m. Furthermore, note that m <fh holds. 
Formula ( fT0| ) represents a discrete time counterpart of 
Formula d7), cf. [20] for details. In particular, a& and 
fifc are monotonically increasing in k, cf. [20, Proposition 
3.3]. Additionally otk and 6tk converge to ax.s and a T <•, 
respectively, cf. [20, Theorem 3.2]. Next, we can use the 
fact that the discrete time suboptimality bounds a.k(N, •), 
fc G No, are monotonically increasing in the parameter m 
and, thus, a^ = ak(N,m) < ak(N,fh) — a,k holds for 
each k G No, cf. [5, Proposition 7.3]. Hence, this inequality 
also holds for the limit, i.e. 

clt,& — lim ak(N,m) < lim a.k{N,rfi) — a T ~ s 

k— >-oo k— foo 

which shows ((9} for T, 5, 8 € Q. ■ 

Note that ay^ defined by Formula |7]) is differentiable on 
(0,T). From Proposition p we can thus conclude a J s " > 
for all 5 € (0, T/2] and ^^ < on [T/2, T). 

Remark 6: Note that the discretization procedure used in 
the proof of Proposition BJ was employed in a purely the- 
oretical fashion. Checking the assumptions on the (discrete 
time) dynamics from [20] is not needed since we only use 
properties of Formula < fT0] > which are independent of its 
connection to discrete time dynamics. 

Combining Corollary |4] and Proposition BJ yields the fol- 
lowing theorem as a direct consequence. 

Theorem 7: Let Assumption |2| and olt,8 > a > hold 
for a minimal control horizon o~G (0, T). Then, ax.s > 5 
holds for all 6 G [<5, T — <5] . Hence, if the stability condition 
a T,(5 > « holds for a desired guaranteed performance bound 
a, then at least the same performance can be guaranteed for 
all control horizons 8 G [5, T — 6] and, thus, also for time 
varying control horizons (8i)i^ C [5, T — 5}. 

Proof: Combining the properties shown in Corollary HI 
and Proposition|5]immediately implies the claimed inequality 
&t,8 > a T,S > 5- The assertion for time varying control 
horizons can be concluded analogously to the stability theo- 
rem for fixed 8, cf. [16], with the same modifications being 
carried out in [5] to extend the discrete time stability theorem 
from fixed to time varying control horizon (Vr(-) is used as 
a common Lyapunov function). ■ 

Remark 8: It is possible to replace the exponential con- 
trollability Assumption [2] by alternative controllability con- 
ditions. In this case, however, in the continuous time setting 
no closed formulas for ot,s as in Theorem [3] are known and 
thus checking monotonocity may become a difficult task. 



Remark 9 ( Sampled-data systems): If a sampled-data im- 
plementation of system (jTJ is necessary, the methodology 
proposed in [20] can be used to determine the required 
sampling rate such that the continuous time estimate is 
approximated arbitrarily well. This allows to transfer our 
results to a sampled-data framework. 

So far, we showed that using larger control horizons is 
not harmful with respect to the performance. Next, for a 
given prediction horizon, we compute parameter combina- 
tions (C, (£) for which stability is ensured by Theorem [3] in 
dependence of the control horizon 8, cf. Fig. 15] 
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Fig. 3. Parameter pairs (C, a) for which «j> g > is ensured by Theorem 
|3]for T = 1. 



The respective structurally different impact of the over- 
shoot bound C and the decay rate /i from Assumption [2] is 
investigated in the next two propositions. 

Proposition 10: Let a decay rate jjl > and times T > 
8 > be given. Then, if Assumption 13] is satisfied with 
a sufficiently small overshoot bound C > 1 the stability 
condition ax,s > is satisfied and, as a consequence, 
asymptotic stability of the MPC closed-loop is guaranteed. 
Proof: Let C be equal to 1. Then, ax.s from |7]) equals 



&T,S = 1 — 



(e^ s - l)(e^ T -V - 1) 



[ e M*( e M(T-«) _ 1)] [ e n(T-S)( e nS _ !)] 

=: 1 - e'" T > 0. 

Since |7]i is continuous with respect to C, choosing the 
overshoot C > 1 sufficiently close to one ensures a^j > 
and, thus, the assertion. ■ 

Note that the performance estimate is bounded by 1— e _AlT 
independent of the overshoot C. Hence, a sufficiently large 
prediction horizon T may be needed in order to guarantee 
a desired suboptimality bound a G (0, 1). Moreover, the 
obtained bound does not depend on the control horizon 8, 

In the following proposition the roles of C and [i are 
reversed. In contrast to the previous observation that, for 
given decay rate /i and prediction horizon T, the maximal 
achievable performance was bounded, here, an arbitrary 
performance a £ (0, 1) can be obtained by choosing the 
decay rate appropriately without prolonging the prediction 
horizon T or changing the overshoot C. 



Proposition 11: Let an overshoot C > 1, times T > 8 > 
and a desired performance bound a € (0, 1) be given. 
Then, if Assumption [2] is satisfied with a sufficiently large 
decay [i > the condition ar,s > a is satisfied. 

Proof: Let the function / : R >0 -> K be defined 
by f{x) — x 1 ^ . Then, its first derivative is given by the 
monotonically decreasing function 

/'(*) = c-aji-i/c- 

Consequently, for real numbers x > x > 1, the inequality 

f(x) — f(x — 1) > /(x) — f(x — 1) holds which implies 

(eMT_l)VCL ( e ^_l)VC > ( e MT)l/C_ (e M)i/c 
( e ^_l)VC*_ ( e ^(T-tf)_ 1 U/C > ^TnI/CL /^(T-Shl/C 

Using these inequalities and the fact that the nominator 
of the second summand in Formula |7]i is smaller than 
( e M<5)i/C( eM (T-<5))i/C = ( e /.Tji/c leads to the estimate 

( e /»T)l/C 



«T,6 > 1 



1 



(eMT)i/o [( eM (T-«))i/c _ !] [( e ^)i/c _ i] 



(e 



jU (T-5)\l/C 



1 



(e 



y.8\l/C 



1 



Since the right hand side of this expression converges to one 
for /i approaching infinity, a (sufficiently large) decay rate /i 
exists such that a^.s > a is ensured. ■ 



Propositions 10 and 11 provide some insight into the 



structure of the proposed controllability Assumption [2] The 
integral L Ce~ tJ,t dt is the destinctive feature. This quantity 
converges to zero for /i approaching infinity. For given /i, 
however, L e _Alt dt is a lower bound. This explains why, 
independently of the overshoot C, the performance guarantee 
cannot be arbitrarily well. Note that Proposition [TT] cannot 
be obtained in a discrete time setting [5]. 

V. Example 

We illustrate our results by computing the a-values in the 
relaxed Lyapunov Inequality |6} along simulated trajectories 
for various control horizon sequences and comparing them 
with our theoretical findings. We consider a continuously 
stirred tank reactor (CSTR) with energy balance and reaction 
A — > B given by the dynamics 

q(x{ 



x\ 



x 2 



x\) 



V 



koX\e 



-E a /x 2 



q(x f 2 



X2) 



-kgxie' 



-E a /x 2 



+ 



a 



X2 



V p ■ c V ■ p ■ c 

where the states x% and X2 denote the concentration in ^ 
and the temperature in K, and the control u corresponds to 
the temperature of the cooling jacket in K, cf. [9], [10]. 

Here, the aim is to stabilize the equilibrium x* = 
(0.5, 350) and u* — 300. The system is subject to the 
physically motivated state constraints X — [0,1] x [0, oo). 
Additionally, control constraints U = [250, 450] are imposed. 
The control u is considered to be piecewise constant on 
intervals of length 8 = 10~ 2 . The optimal control in each 
MPC step is computed over an optimization horizon T = 0.3 



name 


symbol 


quantity 


unit 


flowrate 


q 


100 


m J /sec 


CSTR volume 


V 


100 


irfi 


pre-exponential factor 


k 


7.2- 10 10 


1/sec 


activ. energy [gas const.] 


E a 


8750 


X/K 


heat of the reaction 


h 


5-10 4 


J/mol 


A/B mixture density 


P 


1000 


kg/m 3 


A/B mixture heat capacity 


c 


0.239 


J /(kg ■ K) 


heat transfer 


a 


5-10 4 


W/K 


feed concentration 


} 
x 2 


1 


mol/m 3 


feed temperature 


350 


K 



TABLE I 

Parameters for the CSTR example. 



(although the CSTR can be practically stabilized using MPC 
with T — 0.02, we use the larger T in our simulations in 
order to have more flexibility in choosing 8). 

We use the initial value x = (0.35, 370) and define the 
running costs via 



£(x, u) 



(x - x\f + {x- x* 2 ) 2 + 10- 3 {u - u*f 



which render the state components equally weighted. The 
error tolerance for the used optimization and integration 
method were set to 10~ 7 and 10 -6 , respectively. Moreover, 
we added a truncation region of the stage costs £ of e = 
10~ 12 to compensate for a possible practical stability region 
and numerical errors, cf. [3, Theorem 21] for details. The 
performance index ar.s can then be evaluated along the 
simulated closed-loop trajectory via 

oit,8 = inf otrAn) (11) 

ne{n\3keN :n=kS} 

where the local performance index aT.s( n ) is given by 
Vr(xn(n5;xo)) - Vr(x fl ((n+ 1)6; xq)) 



a T ,s(n) = 

j{l{x^,{t; x^{n8; x )),fx(t; x^nS; x ))))dt - e 
o 

with \i = ht.s if the denominator of the right hand side is 
strictly positive and aT,s( n ) — 1 otherwise. Note that ar,s 
is negative if the value function increases. 

Fig. [4] shows the closed-loop solution in dependency 
of different fixed and time varying control horizon se- 
quences. Although the trajectories corresponding to 8 G 
{0.15, 0.25, 0.30} appear to be less efficient due to avoiding 
the turning point close to (0.4, 325), a computation of the 
performance index ao.3,<5 reveals the opposite, cf. Fig. B] If 
the control horizon is increased, the performance index rises 
as expected from our theoretical results in Section llV] 

For this particular example, symmetry of ar.s with respect 
to the control horizon 8 is not observed. Instead, the per- 
formance index increases (almost) monotonically. There are, 
however, examples in which the suboptimality degree decays 
rapidly for longer control horizons, cf. [7]. Utilizing the out- 
come of Fig. Bl we restricted 8 to [0.1, 0.3] and computed the 
closed loop solutions shown in Fig. HI From this, we observed 
that the suboptimality bound arj. 3,0.1 ~ 0.3346 (obtained for 
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Fig. 4. Closed-loop solutions for different fixed control horizons <5 and 
400 closed-loop solutions for randomly time varying <5 £ [0.1,0.3]. 
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Fig. 5. Development of ciq.3 S f° r fixed control horizons 5. 00.3,0.1 yields 
a lower bound for time varying control horizon S 6 [0.1, 0.3]. 



fixed control horizon 6 = 0.1) holds for control horizons 
randomly varying in a rather large interval. For short control 
horizons on the other hand, the relaxed Lyapunov Inequality 
(|6} — our main tool in order to ensure stability — is violated, 
cf. Fig. 15] In an a posteriori analysis, we found that the 
closed-loop costs J^ c (xq) are almost constant with respect 
to S. In conclusion, this example demonstrates that packet 
dropouts and non-negligible delays can be compensated via 
time varying control horizons without affecting our (nominal) 
stability and performance estimates. 

VI. Conclusions 

We presented stability and performance estimates for 
nonlinear model predictive control with time varying con- 
trol horizon in a continuous time setting. In particular, 
we deduced symmetry and monotonicity properties for the 
stability condition introduced in [16] which were exploited 
in order to show that no additional assumptions are needed 
compared to the analysis of schemes with fixed (short) 
control horizons. The results can be used in order to obtain 



(nominal) stability and performance estimates for prediction 
consistent networked control schemes. 
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